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Abstract 

In this paper we prove the continuity of stable subspaces associated 
to parabolic- hyperbolic boundary value problems, for limiting values 
of parameters. The analysis is based on the construction performed in 
|MZ| of Kreiss' type symmetrizers. 



1 Introduction 

This paper is motivated by the analysis of multidimensional small viscosity 
boundary layers or shock profiles (see |MZ], jGMWZlj |GMWZ2| ). The 
stability criterium for layers or profiles is described through the nonvanishing 
of an Evans function. After freezing the tangential coefficients, performing 
a tangential Fourier-Laplace transform and using a reduction to a normal 
form, the problem is reduced to an ordinary N x N constant coefficient 
differential system in the normal variable: 

(1.1) d x U = G{p,()U + F for x>0, 

where p are parameters and £ denote the tangential frequency variables 
(see (HZ1)- Typically, ( = (r - i~f,v) G C x R^ 1 where t — i-y is the 
Fourier-Laplace frequency dual to time t and r\ = (771 , . . . , rjd-i) are the 
Fourier frequencies dual to the spatial tangential variables y. The equation 
is supplemented by boundary conditions 

(i-2) r{p, qu {x=0 = . 
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For p = |e| > and 7 > 0, G(p,() nas no purely imaginary eigenvalues 
and the bounded solutions of the homogeneous equation (|1.1|) are 

U(x) = e xG{p '°U(0) , U(0) G E_(p, C) 

where E_(p, £) is the stable subspace at +00 of Q{p, (), that is the space 
generated by the generalized eigenvectors associated to eigenvalues belong- 
ing the the half space {Re^ < 0}. The stability condition, that is the well 
posedness of the boundary value problem <|l-2j) in L p spaces, reads 

(1.3) C 7V = E_(p,C)0kerr(p,C), 
or 

(1.4) D(p,() :=det(E_(p,C),kerr(p,0) #0 

where, for linear subspaces E and F of C^, det(.E, F) is zero if dmiE + 
dimi* 1 7^ A?" and equal to the determinant obtained by taking orthonormal 
basis in E and F when dimE + dimF = N; this is necessary for mul- 
tidimensional stability. The function D is commonly known as an Evans 
function. 

The uniform stability condition requires that D (or a rescaled version for 
high frequencies) be uniformly bounded from below by a positive constant 
when p > and 7 > 0. This has been shown in several contexts to be 
sufficient for multidimensional stability; see [Z], [MZ],[GMWZl-2]. 

Of particular importance is the analysis for low frequencies where the 
hyperbolic character becomes crucial. There one uses polar coordinates 

(1-5) ( = p(, P = \C\- 

The main result of this paper is: 

Theorem 1.1. Under Assumvtion \3.1\ below, the linear bundle E(p, p) := 
K—(p,pQ has a continuous extension to p = 0, 7 > 0. 

This result was known under additional assumptions such as strict hy- 
perbolicity of the hyperbolic part or constant rank of branching points, see 
|Rouj . |ZSj . or |Z]. The Assumptions 3.1 are the somewhat minimal condi- 
tions of hyperbolicity and parabolicity as in [MZ . An important corollary 
is that 

(1.6) D(p,p() = /3A(p,() + o(l), as p^0 
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uniformly for \Q = 1 with 7 > 0, where A(p,Q is the Lopatinski determi- 
nant of the appropriate limiting hyperbolic boundary value problem, and (3 
is a transversality coefficient for the inner layer ordinary differential equation 
(the viscous shock or boundary profile problem); indeed, the two statements 
are essentially equivalent. This reduces the somewhat abstract uniform sta- 
bility condition, at least in the crucial low frequency regime, to a pair of 
checkable conditions, associated respectively with outer (hyperbolic) and 
inner (stationary parabolic) problems in the formal matched asymptotic ap- 
proximation of the boundary layer, and is one of the central conclusions of 
[ZS],[Z],[Rou]. 

In the original planar shock stability analysis of [Z], which was carried 
out by explicit estimates on the Laplace-Fourier inversion formula, con- 
tinuous extension of stable subspaces was also used in an essential way in 
obtaining the linearized decay estimates from which nonlinear stability ulti- 
mately follows. In the more recent curved boundary layer analysis of [MZ], 
carried out by Kreiss symmetrizer techniques, it was shown that the uniform 
stability condition alone suffices, whether or not continuity holds, yielding 
applications to much more general situations than considered in [Z] ; this was 
possible because of the flexibility of the symmetrizer method. This abstract 
result, however, sidestepped the issue of checkability of the uniform stability 
condition/validity of the fundamental relation ()1.6|) . 

With the establishment of Theorem we now resolve this issue com- 
pletely, under the same minimal hypotheses considered in [MZ] . The matrix 
perturbation problem associated with determination of stable subspaces is 
a rather complicated one involving many parameters, and would be difficult 
if not impossible to describe completely in this generality. In the inviscid 
case, this problem can be sidestepped thanks to the block structure condi- 
tion enjoyed by systems with constant multiplicity characteristics; see [Me3]. 
However, the additional parameter p generically destroys this block struc- 
ture, and so the issue apparently arises again. A notable aspect of [MZ] 
was the construction of Kreiss symmetrizers without the use of block struc- 
ture in either 7 or p: that is, the identification of different roles played by 
Laplace and Fourier frequencies in the analysis. Here, we show that, some- 
what surprisingly, the same symmetrizer construction used in [MZ] to avoid 
the need for continuity of subspaces in fact establishes that continuity! That 
is, the structure of stable subspaces is already hidden in the details of the 
nonstandard symmetrizer construction. 

As the proof of the main theorem relies on the use of symmetrizers 
for Q(piOi we fi rs t review the construction performed in |MZj of Kreiss' 
type symmetrizers for parabolic-hyperbolic boundary value problems and 
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we prove the new continuity theorem next. 



2 Symmetrizers 

Recall first the essence of the "method of symmetrizers" as it applies to 
general boundary value problems on the half line {x > 0}: 

(2.1) d x u = G(x)u + f, r«(0)=0. 

Here, u and / are functions on [0, oof values in some Hilbert space Ti, and 
G(x) is a C 1 family of (possibly unbounded) operators defined on T>, dense 
subspace of Ti. 

A symmetrizer is a family of C 1 functions x i— > S(x) with values in the 
space of operators in Ti such that there are Co, A > 0, 5 > and C\ such 
that 

(2.2) Vx, S(x) = S(x)* and \S(x)\ < C , 

(2.3) Vx , 2Re S(x)G(x) + d x S(x) > 2AId , 

(2.4) S(0) > 6Id - cir*r . 

In (|2.2|) . the norm of S(x) is the norm in the space of bounded operators in 
Ti. Similarly S*(x) is the adjoint operator of S(x). The notation ReT = 
|(T + T*) is used in (|2,3[) for the real part of an operator T. When T 
is unbounded, the meaning of Re T > A, is that all u £ T> belongs to the 
domain of T and satisfies 

(2.5) Re (Tu,u) > X\u\ 2 . 

The property (|2.3j) has to be understood in this sense. 

Lemma 2.1. // there is a symmetrizer S, then for all u € Cq([0, oo[;TL) D 
C°([0,oo[;£>) ; one has 

r 2 

(2.6) A||«|| 2 + %(0)| 2 < ^\\d x u - Gu\\ 2 + Ci|ru(0)| 2 . 

A 

Here, | • | is the norm in Ti and || • || the norm in L 2 ([0, oo[; Ti). 

Proof. Taking the scalar product of Su with the equation 1)2.1(1 and inte- 
grating over [0, oof, (|2.2|) implies 



-(S(0)u(0),u(0)) = / d x (Su,u)dx 
(2-7) r J 

((2ReSG + d x S)u,u)dx + 2Re / (Sf,u)dx. 
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By Q, 

((2Re SG + d x S)u,u)dx > 2A||u|' 2 
By (|2.4|) and the boundary condition Tu(0) = 0, 

(5(0)«(0),«(0)) > 5\u(0)\ 2 - Ci|rn(0)| 2 . 

By d2U) 

(5/,«)dx| < 2Cb||/|| HI < ^||/|| 2 + A|M| 2 . 

Thus the identity (|2.7|) implies the energy estimate (|2.fij) . □ 

In applications to evolution boundary value problems, is a tangen- 
tial differential system G(t, y, x, dt,d y ). One usually studies (|2.1|) in weighted 
spaces e 7 *L 2 , with 7 > for forward propagation. With u y = e _7 *u and 
/ 7 = e~ 7 */, the equation (|2.1j) is transformed into 

(2.8) d x u^ = G 7 (x)n 7 + / 7 , Tn(0) = . 

with G 7 (t, y, x, dt,d y ) = G(t, y, x, dt + 7, d y ), which is studied in I? . When 
G has constant coefficients in (t,y), the equation (j2.8|) transformed into 
a system of o.d.e. in x using a tangential space time Fourier transform. 
For "slowly" varying coefficients, one first freeze the coefficients and then 
perform the Fourier transform, yielding a family of systems depending on 
the parameters (t,z) and the frequencies C = C^^t): 

(2.9) d x u.y = G(t, y, x, £)n 7 + / 7 = , ru 7 = , 

where G(t, y, x, £) = G(t, y, x, it + 7, irj). 

In this framework, one can look for symmetrizers for the system (|2.9j) . 
This is a matrix S(t, y, x, () such that 

S(t,y,x,() = S(t,y,x,() and \S(t, y, x, ()\ < C , 

(2.10) V 7 >0, 2ReS(t,y,x,0G(t,y,x()+dJ(t,y,x,() > jld, 
5(t,y,0,C) > Sld-CtF*?. 

When G and S do not depend on (t,y), one immediately gets symmetrizers 
for (|2,8[) . taking S as Fourier multiplier. When S depends on (t,y), one can 
try to take as a symmetrizer the pseudo-differential operator 

5 7 (x) = S(t,y,x,D t ,D y ,j) 
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or a para-differential version of it. The problem is then to convert the sym- 
bolic estimates ()2.10(l into estimates for the operators S(x). This strategy 
has been shown to be successful for hyperbolic problems ( |Kr j . |CPj ^ and 
small viscosity parabolic-hyperbolic problems jGMWZlj . |GMWZ2j ). 

A further reduction to constant coefficients in x can be performed when 
G converges at an exponential rate to a limit when x tend to infinity: 

(2.11) \G(x,Q-G(oo,C)\ <Ce~ dx 
for some 6 > and £ in a neighborhood of £. Then: 

Lemma 2.2 (|MZJ). There is a neighborhood uo of Q and there is a matrix 
W defined and C°° on [0, oo[xu; such that 

i) W _1 is uniformly bounded and there is 0\ > such that 

(2.12) |W(x,C) -Id| < Ce- dlX ) 

ii) W satisfies 

(2.13) 9 x W(x, C) = G(x, C)W(x, C) - C)^(oo, C) • 
The substitution 

v {x) = W-\x, ()u 7 (x) , g(x) = W-\x, (fyx) 
transforms l)2.9j) into 

(2.14) d x v = G(oo,()v + 9, ri(CMO)=o, 

with Ti(c) := rw-^o.C)- 

All this motivates the construction of symmetrizers for constant coeffi- 
cients systems (|l.l(l with boundary conditions (|1.2j) . Note that the reduction 
to constant coefficients above leads to frequency dependent boundary condi- 
tions. In the sequel, we concentrate on hyperbolic-parabolic systems arising 
in the study of boundary layers or shock profiles. 

3 The main result 

Consider a first order N x N system 

d 

(3.1) L 1 {jp,d)u:=d t + Y, A M d 3 
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The variables are t G M., y G W d ~ l and x G R+ :=]0, +oo[. Moreover, 
<9j = <9 yj . for j G {1, . . . , d — 1} and = d x . In addition, p G M M is a 
place holder for the variables themselves, the unknown in case of quasilinear 
systems and all the other parameters such as source terms or control. 
Next, we consider a parabolic viscous perturbation of (|3.1|) 

(3.2) L(p,d)=L 1 (p,d)-L 2 (p,d) 
with 

(3.3) L 2 (p,d):= d^Bj^dku) . 

l<j,k<d 

This kind of systems arises in the study of boundary layers or shock profiles, 
as linearized form of viscous perturbations of hyperbolic conservation laws, 
see e.g. jMZ] I^MWZlj . |GMWZ2| . In [MZj |GMWZ2| we considered rather 
the small viscosity problem L\ — eL 2 . Rescaling the variables as t = et, y = 
ey,x = ex reduces to s = 1 and the rescaled variables are much more 
adapted to the geometrical analysis of the symbol. 

Assumption 3.1. 

(HO) The Aj and Bj k are C°° functions of p on the open subset O of 
M.)M, with values in the space of N x N real matrices, 

(HI) For all p G O, the eigenvalues of ^2^jAj(p) are real and semi- 
simple and have constant multiplicities for p G O and ^ G M. d \ {0}. 

(H2) There is c > such that for all p G O and £ G M. d the eigenvalues 
°f i Ej=i £j A ji u ) + Ej,fc=i £j€kBj,k(u) satisfy Re/i > c|£| 2 . 
(H3) For allpe O, there holds det A d {p) ^ 0. 

We perform a Laplace-Fourier transform with respect to t and a Fourier 
transform with respect to the variables y. This amounts to performing the 
substitution 

(3.4) d t i-> ir + 7 , d y . h-> ir]j 

in the definition of L, with r £ R, 7 £ M and rj = (771, . . . ,i]d-i) £ 
In applications, to study the forward propagation, dictated by the sign of 
L 2 , one restricts 7 to be nonnegative. After this Laplace- Fourier transform, 
L becomes a second order differential operator in x depending on the pa- 
rameters p and C = ( r ; r ?)7) € We reduce this system to first order 
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introducing the variable U = (u, d x u) with values in M. 2N . This leads to the 
following 2N x 2N system 

(3.5) C:=Idd x -g(p,(), G(p,0=(m ! j 
where 

/ d-l 
A(p, C) = {B d4 )~ l [A d - iVj(B jid + Bdj)) 

i=i 

d-l d-l 

M(p, () = {B d4 )- 1 [{ir + 7) + ir lj A j + Yl VjVkBjjt 

j=l j,k=l 

and the matrices Aj and Bj^ are evaluated at p. 

The Assumption 13, II immediately implies the following result (see |MZj 
for a detailed proof). 

Lemma 3.2 (Spectral analysis of Q). 

i) When ( 7^ and 7 > 0, 9(p, C) Af eigenvalues, counted with their 
multiplicities, in Re /i > and iV eigenvalues in Re < 0. 

ii) When C = 0, G(p, 0) has as a semi-simple eigenvalue, of multiplicity 
N. The nonvanishing eigenvalues are those of (-B^d) -1 A d . 

For ( / and 7 > we denote by E_(p, £) [resp. E + (p, £)] the spec- 
tral subspace of <?(p, C) associated to the spectrum in the left hand half 
space {Re fi < 0} [resp. {Re fi > 0}]. Both have dimension ./V and depend 
smoothly on p and £ as long as £ 7^ and 7 > 0. To analyze the behavior 
of this spaces when £ tends to zero, we introduce polar coordinates C = PC 
with /o = |C| and set 

(3.6) E(p,C,p):=E_(p,pC), forp>0. 
We can now state the main result of this paper: 

Theorem 3.3. Under Assumption lff.il the linear bundle E(p, £, p) /ias a 
continuous extension to p = 0, £ € \ {0} rnf/i 7 > and p £ O. 

Remark 3.4. As already said, this result was already known under addi- 
tional hypothesis |Rouj . |ZSj . [Z]. In |MZj we used that the E_(p, 7, p) is 
continuous up to p = when 7 > and that E(p, £, 0) has the form 

(3.7) E(p, C, 0) = E^ p (p, C) E p _ ar (p, 0) , 
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see section 4 below. Moreover, E_ yp is the negative subspace corresponding 
to the hyperbolic system L\. Following Kreiss |Krj (see also j(JPj and |Me3j 
for the extension to hyperbolic systems of constant multiplicity) Kj^ p (p, () 
can be extended continuously up to 7 = 0. Therefore, in |MZj we used the 
existence of the double limit 

lim lim E_(p, (, p) 

to establish desired linear decay estimates. The Theorem 13.31 asserts the 
existence of the limit as 7 and p tend to zero, independently. 

Remark 3.5. Consider a boundary condition T(p, QU = for ()3.5|) where 
r is a smooth function of p and £. The weak stability condition is given by 
(|1.4|) and the strong stability condition reads 

(3.8) \D(p,0\ = |det(E_(p,C),kerr(p,C))| >c 

for some constant c > 0. When it holds, passing to the limit in p implies 
that I det(E_(p, 0), ker T{p, 0))| > c for 7 > 0. Next, passing to the limit 
in 7 shows that this condition holds up to 7 = 0. Thus, 

(3.9) C 2N = E_(p,C,0) 0kerr(p,O) 

when 7 = 0. Theorem 13.31 when 7 = and the known smoothness when 
7 > 0, imply that conversely, if the transversality 13.91 holds at some point 
(p,C) with C 7^ an d 7^0, then the uniform stability condition 13.81 holds 
for (p, £, p) in a neighborhood of (p, (, 0) with 7 > and p > 0. 



4 Strategy of the proof 

By Lemma 13.21 for £ small the eigenvalues of Q split into a group of N 
eigenvalues close to and a group of N eigenvalues away from the imaginary 
axis. More precisely, given a point p, there is a C°° invertible matrix V 
defined for (p, £) in neighborhood uj of (p, 0) such that Q\ := V~ l QV has the 
block diagonal form 

H(p,C) 

p( P ,o 

with H(p,0) = and P(p,0) = (-B^) -1 ^- The eigenvalues of P(p,C) 
satisfy |Re /i| > c for some c > and 

d-i 

(4.2) IT = -(Ajr 1 (ir + 7)Id + X;^i) + 0(|C| 2 )- 

i=i 



(4.1) ai(p,C) 
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Therefore, for (p, () € uj, with C 7^ an d 7 > 0: 

(4.3) E_(p,C) = V(p J C)(E?(p ) C)®l-(p,0) 

where [resp. Ef] is the spectral space associated to eigenvalues of H 
[resp. P] in {Re fi < 0}. 

On one hand, because the spectrum of P does not intersect the imaginary 
axis, the definition of the spaces E^ extends smoothly to (p, Q € uj. On the 
other hand, the eigenvalues of H are small, and we now investigate how they 
split in the complex complex domain when £ is small. In polar coordinates, 
there holds 

(4.4) H(p, pC) = P H(p, C, p) = p{H (p, C) + 0{p)) 
with 

d-l 

(4.5) H (p, C) = -{Ad)' 1 (if + 7 )Id + Y, ivAj) 

3=1 

Therefore 

(4.6) E_(p,pC) =V( P ,p()(E / [(p,(,p)®eF(p,p()) 

where E^ is the negative space associated to H. 

The hyperbolicity assumption (HI) implies that Ho(p,Q has no eigen- 
values on the imaginary axis when 7^0. This remains true for H(p,£, p) 
for p small enough (depending on 7). In particular, for all £ £ M. d+l with 

7 > 0, the definition of the bundle K H (p, £, p) extends smoothly to a neigh- 
borhood of (p, c, 0) in R M x R d+1 x R. This remains true for E(p, p£) since 
E^(p, pC) an d V(p, p() are smooth in p£. 

Thus, it remains to prove Theorem 13.31 on a neighborhood of a point 
(p, C, 0) where 7 = and C / 0. 

The strategy of the proof is as follows. We first construct explicitly 
spaces E_ and E + such that 

(4.7) C 2N = E_ E + , dimE_ = dimE + = N . 

Indeed, extending Kreiss' analysis to hyperbolic systems of constant multi- 
plicity as in |Me3j . one can show that the bundle E^°(p, £) extends contin- 
uously to 7 = and we take E_ to be equal to the value of this extension 
a ^ (P)C)> but we will not use this result, see below. Next the main point 
is to construct symmetrizers. In the next statement we denote by II_|_ the 
projectors on E ± associated to the decomposition (|4.7j) . 
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Theorem 4.1. For all (p,Q with ( 7^ and 7 = 0, there are spaces E ± 
satisfying 1)4.71) and such that for all K > 1 there are a neighborhood Q of 
(p, C) in R M x R d+1 x R, a C°° mapping S from O to i/ie space 0/ 2iV x 2iV 
matrices and a constant c > such that for all (p, p) £ f2, 



(4.8) S(p,£,p)=S*(p,C,p), 

(4.9) (S(p,(,p)U,U) > k 2 \U + U\ 2 - \n_u\ 

and for all (p, p) £ fi wii/i /? > and 7 > 0: 

(4.10) Re(S(p,C,p)0(p,pC)tf,CO >cp(7 + p)|f/ 



Proof of Theorem \3.3\ assuming Theorem 14.11 

a) Consider k > 2 and f2 given by Theorem (|4. 1() . For (p, p) G f2 with 
7 > and p > 0, let U G E_(p, /<)• Then 

is exponentially decaying at +00 and satisfies d x ll = G(p, pC)U ■ Therefore 
multiplying by S(p,C,p) and integrating by parts as in Lemma 12.11 yields 
thanks to (jHHJ): 

(SU,U) + 2Re / (SGU(x),U(x))dx = 
Jo 

By (|4.10|) . the integral is nonnegative. Therefore (SU,U) is nonpositive 
which by l|4~5J) implies that k\U + U\ < \U_U\. Thus 

(4.11) V[/GE_(p,pO : \U.U\<—^—\U\. 

K — 1 

This implies that the mapping n„ from E_(p, into E_ is one to one 
and since both spaces have dimension N, it is a bijection. Therefore, there 
is a mapping <3?(p, (, p) from E_ to E + such that 

(4.12) E-0,/<) = {u + A(p,(,p)u :u£E_} 
and 

(4.13) VueE_: \A(p,(,p)u\ < Id. 

k — 2 

Since k is arbitrarily large, this proves that 

(4.14) E_ = lknE_(p,pC), 

where lim means that (p, (, p) tends to (p, C, 0) with p > 0, 7 > 0. 
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b) Since this is true for all (p, Q this proves that E_ (p, f , p) := E_ (p, pQ 
extends to points (p, C, 0) with 7 > 0. We denote by E_ (p, (, p) this extension 
and (|4.14|) implies that for all (p, £) with 7 > 0, there holds 

(4.15) E_(p,C,0) = limE_(p,C,p). 

Consider again a given point (p, C) with 7 = 0. For k > 2, let be 
given by Theorem 14.11 For all (p',C',0) ^> thanks to (|4.15|) . we can let 
(Pi Ci P) tend (p', v , 0) in the sense of lim. Therefore, (|4.11l) implies that for 
all (p, C, p) G ^ with 7 > and p > 0: 

V[/GE_(p,C,p) : |n + c/|<^-|c/|. 

K — 1 

Arguing as before, this implies that 

E_(p,C,0) =E_ =limE_(p,C,p), 

where the limit is taken for (p, p) tending (p, 0) with p > 0, 7 > 0. 

This means that the bundle E_(p, p) is continuous in (p, p) for 7 > 
and p > 0. □ 



5 Proof of Theorem 14.11 

Though not stated explicitly, the result of Theorem 14. II is given in Appendix 
A, section 3, of MZ . For the convenience of the reader, we give here the 
main steps of the proof. 

a) We first remark that it is sufficient to prove the theorem for k large 
enough. Next, we can replace (|4,9[) by the weaker condition 

(5.1) (SU, U) > k 2 \U + U\ 2 - \U_U\ 2 , 

for 5 = S(p, C,0). Indeed, if (|5.1|) holds, then, by continuity of S, on a 
possibly smaller neighborhood of (p, 0), there holds 

(s(p,c,p)u,u) > \k 2 \U + u\ 2 - 2\n_u\ 2 , 

and therefore satisfies (|4,8|) for k/2 and also (|4.8|) and (|4.1()j) with the 
constant c/2. 
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b) Recall that Q\ = V l QV has the form l|4.1jl . We construct spaces 
E± and E^ such that 

(5.2) C N = E H , C N = E P ®E P dimE H + diniE^ = N. 

Denote by ITjf and the corresponding projectors. Suppose that for all k 
large enough there are symmetrizers S H and S p defined on some neighbor- 
hood of (p, () and such that 

(5.3) S H = (S H )* 

(5.4) (s H u, u)) > k 2 \u p u\ 2 - \u p u\ 2 , 

(5.5) Re S H H > 0(7 + p)Id , for p > , 7 > , 
and 

(5.6) S p = (s p y 

(5.7) {S P U, CO) > K 2 \U P U\ 2 - \u p u\ 2 , 

(5.8) Re S P P > eld . 

Here, we have denoted by S_ H and S_ p the value of S H and S p respectively 
at the base point (p, £, 0). 

Then we choose E ± = V(Ef §E^) and S = (V^ 1 )*^"^" 1 with 

(5-9) Si(p,C) = f( 5 * 5 °p) 

By and (|53|) . S is clearly self adjoint and (f53|) and (|S~g)) imply (|4~TU|) . 
In addition, for U = V{U H © U p ) one has 

(St/, 17) = 5{S H U H , U H ) + 5(S P U P , U p ) . 

Because n ± = V(n± ©n^)V _1 , there is C such that 

|n^| 2 + |n^c/ p | 2 >^|n + c/| 2 
|n_c/| 2 <c(|n^| 2 + |n p c/ p | 2 ) 

Therefore, choosing 5 = 1/C, (j5.4|) and (j5.7|) imply (|5.1j) with k replaced by 
k/C and it is sufficient to construct the symmetrizers S and S p . 
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c) The construction of S p is easy. Because the eigenvalues of P do 
not belong to the imaginary axis, one can perform a further reduction (or 
choose V in (|4.1|) 1 such that for (p,() m a neighborhood of (p, 0): 



p(p,0 



p+(p,C) o 
o P-(p,C) 



where the spectrum of P± is contained in {±Re/x > 0}. In this basis, the 
spaces ^±(p, C) are constant, we call them E ± . They are generated by the 
first and last vectors of the basis corresponding to the blocks P+ and P_ 
respectively. 

There are self- adjoint matrices S± and c > such that 
±ReS£P± > eld, S£>Id, 5 p <Id. 
For instance, one can choose with appropriate positive constants C± 



si = c± 

Then the symmetrizers 

s p ( P ,C) = 



TtP±\*TtP± 



dt. 



i 2 sP(p,() o 
o -s*(p,0 

satisfy ([Ojl to (JEHJ). 

Remark 5.1. The symmetrizers S constructed above are smooth func- 



tions of (p, (). This was important in |MZj in the quantization of the sym- 
bolic calculus. 

d) We now turn to the construction of S H . As already mentioned 
in (|4.4|) . H is a perturbation of Hq. We denote by fj, = the distinct 
eigenvalues of Hq := Hq(p, Q located on the imaginary axis. There is 5 > 0, 
a neighborhood f^o of (p, C, 0) and a smooth matrix V(p, C, p) such that 



(5.10) 



V~ lr HV 







Hi 



such that each block 7ik has its spectrum either in {|Re/i| > 25} or in the 
ball of radius 5 centered at Moreover, we can assume that the balls of 
radius 25 centered at the do not intersect each other. 
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We denote by Nk the dimension of the block 7ik- By Assumption 13,11 
the 7ik have no eigenvalues on the imaginary axis when p > 0, 7 > and 
p + 7 > 0. Therefore, the number of eigenvalues of Hk in {±Re/i > 0} is 
constant for p > 0, 7 > and p + 7 > 0. We denote it by iVj^-t. Because 
the total number of eigenvalues of Q in {Re p < 0} is N, we have 

(5.11) J2 N k- + dimE p = N. 

Arguing as in step b), it is sufficient to construct spaces E± satisfying 

(5.12) C Nk = E* © , dimE fc = N k ,- , 

and for all k large enough, symmetrizers S k , C°° on a neighborhood of 
(p, C, 0), such that 

(5.13) S k = (S k )* 

(5.14) (S k U, U)) > K 2 \U k + U\ 2 - \U k _U\ 2 , 

(5.15) Re S k H k > 0(7 + p)Id , for p > , 7 > . 

In (|5.14j) . S k = Sk(p,(,0) and n^. are the projectors associated to the 
decomposition (|5.12l) . 

e) If the spectrum of 7i k lies in {Re^ > 25} [resp. {Re/i < —25}, then 
arguing as in step c), we set E*i = {0} [resp. E^ = C^] and S k = K 2 S k 
[resp. S k = —S k ] where S k {p,C,,p) is a self adjoint matrix, defined and C°° 
on a neighborhood of (p, 0) and such that Re S k > Id and Re S k 7i k is 
positive definite [resp. KeS k < Id and KeS k 7i k is negative definite]. 

f) Consider now the case where the spectrum of TL k is contained in the 
ball of radius 5 centered at p k = i£ fc . We note that i£ fc is an eigenvalue of H_ Q 
if an only if —f = X(p, 17, £ ) where X(p, r/, £) is one among the eigenvalues 

Suppose first that £ fc is a simple root of the eigenvalue equation, that is 

(5.16) dz\(p,v,i k )^0. 

Denote by the multiplicity of A (see Assumption (HI)). Then, according 
to Mr.'! , TL k {p, C, 0) is a scalar matrix q k (p, C)Id and by Lemma 2.10 of MZ 

(5.17) n k ( P X,p) =q k (pX)U + pK k (p,(,p). 
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Moreover, qk is purely imaginary when 7 = 0, qk ■= djReq k (p,() does not 
vanish and qf.KeTZ k (p, 0) is definite positive. 

When <jk > [resp. qk < 0], for 7 > small, the eigenvalue % has a 
positive [resp. negative] real part. Thus we set E_ = {0} [resp. Eji = C Nk ]. 
Next we choose S k = n 2 S k [resp. S k = —S k ] where S k (p,(,p) is a self 
adjoint matrix, defined and C°° on a neighborhood of (p, 0) and such that 
ReS* > Id and ReS k K k is positive definite [resp. Re5 fe < Id and ReS k TZ k 
is negative definite]. 

g) We now come to the most difficult part of the construction, when £ fc 
is a multiple root of the eigenvalue equation, that is, when there is > 2 
such that 



%A = ...^ fc_i A = 0, ap^o 

a t (P; ?7> (Note that, since A is real analytic in £, there is always such an 
integer The case ^ = 1 is 1)5. 16|) ). 

^From |Krj in the strictly hyperbolic case and from |Me3j for the exten- 
sion to constant multiplicity, we know that there is a smooth matrix Vg (p, C) 
such that 



(5.18) Q k (p,C) := (V k )- l (pX)H k (pX,0)V k ( P ,C) 



Qk 











where the subblock Qk is v\~ x V}. matrices and there are such blocks, 
where is the multiplicity of the eigenvalue A. In particular the dimension 
of TL k is Nk = ctkVk- Moreover, 



(5.19) 



Q k -=Q k (p,C) 



2* 




1 








1 



In addition, with Ralston's lemma |Ralj . we can assume that only the first 
column of Qk does not vanish: 



(5.20) 



Qk(p,(,o) 



qk 



0...0 

0...0 
0. . .0 



where qk denotes the lower left hand corner of Qk, that Qk has purely 
imaginary coefficients when 7 = and that qk := d^Reqk(p, C) / 0. 
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Next from Lemma 2.10 of [MX; . we know that there is a smooth matrix 
V k on a neighborhood of (p, 0), which extends Vo, such that 

(5.21) U\ := (V k y l H k V k = Q k (P,0+pK k (P,Cp), 

where the matri: 
of size Vf. x v\~ : 

(5.22) 



where the matrix 1Z k has the following decomposition in x «(. blocks R k q 



T>k 
K 1A 



r>A; 

K l,a k 



Moreover, only the first column of the R k q does not vanish: 



(5.23) 



0...0 

0...0 
0...0 



where r k q denotes the lower left hand entry of R k ^ q - In addition, denoting 
by R\ the x matrix with entries r k q , the matrix ^Re R\(p, (, 0) is 
definite positive. 

Arguing as in step a) and b), we are reduced to construct spaces and 
symmetrizers for each block Ti\. 

h) Prom |Krj |CPj . or simply from a direct analysis of the model case 



(5.24) 



— A- 

o i 



i 



—A' 







* 


o . 


. o " 







* 


o . 


. 


+ 7 








1 






o . 


. 


ik. 




. qk 


o . 


• . 



we know that the candidate for the limit negative space for this block is the 
space generated by the first (5k vectors of the canonical basis : 



(5.25) 
where 

(5.26) 



E*(p,C) = C A x {0}"*- 



0k 



when Vk is even . 



\(yk ± 1) when Vk is odd and =F Qk > . 
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We also introduce 



Thus 



E^(p,C) = {0}^ 



y*k = E_ © E* 



In the block decomposition of C k into a k factors C fc , let 
E£(p,C,0) =E^(p,C)©---©Ei(p,C). 

Note that is exactly the number of eigenvalues in {Re p < 0} of the model 
Q5.24J1 for 7 > 0. This is still true for Qkip, £> ?7>7)- Thus the number of 
eigenvalues in {Re /i < 0} of Q fc is Vkh when 7 > and therefore the number 
of eigenvalues of TL\ in {Re p < 0} when p > 0, 7 > and p + 7 > 0, which 
is constant on a neighborhood of (p, 0), is 

N k = v k f3 k = dimE fe . 

Therefore, the condition ()5.12|) is satisfied. 

We now proceed to the construction of the symmetrizers. We construct 
S k in the block decomposition of 7i\ 



(5.27) 



with 



(5.28) 



S k 
S k 



S k {p, C, p) = E k + E k (p, C) - i"fF k - ipF' k , 



where E k and E k are real symmetric matrices, and F k and H k are real and 
skew symmetric. E k is constant and has the special form 









e k ,i e k2 







e k ,i 
e k ,2 



ek. 



"k 



Moreover Ek(p, C) = and Fk and F' k are constant. 

The order of the construction is as follows. One first chooses E k , Ek and 
F k as in |Krj (see also |CPj ) to construct symmetrizers for Q kl that is for 
p = 0. The new part lies in the choice of F' k . 
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1. Choice of E k . First choose the real coefficient e k) \ such that 
(5.29) e M g fc > 3. 

Next, the coefficients e k i for I > 1 are chosen successively to achieve that 
there is c > such that 



(5.30) (E k w,w) > c(k 2 \7tXw 



\-K k _w\ 2 



where tt_± is the projection onto E± in the decomposition C Uk = E^ © E_. 
(cf Lemma 5.6 in |CPj . Chap. 7). 

Next, following jCPj (cf (equation (5.5.3) in Chap. 7), (j5.29|) implies 
that there is a constant C and a neighborhood of (p, () such that 



12 nw'\\ 



(5.31) Re(E k d 1 Q k (p,Qw,w) >2|wi| 2 -C| 

with w\ the first component of w € <£ Uk and w' € C^ fe_1 denotes the other 
components. 

2. Choice of E k . Recall that Q k = i(£ Id+iV^) where N k is the Jordan 
matrix of size v k (see (j5.19|0 . The form of E k is chosen so that E k (^ k Id+N k ) 

is real and symmetric . Next, the real matrix E k (p, Q is chosen so that such 
that (E k + E k )(iQ k ) is real and symmetric when 7 = 0. This is achieved in 
|Krj |CP| using the implicit function theorem and the property that \Q k is 
real when 7 = 0. 



3. Choice of F k . Following Kr CP , there is F k real and skew sym- 
metric such that 

Re (F k N k w, w) > -\wi\ 2 + (C + l)\w'\ 2 . 
where C is the constant in (|5.31jl . As a consequence, we have 
Re ((E k + E k - i^F k )Q k ) = jD k , 

(5.32) 

D k (p,C) = Re (E k d~,Q k (p, + Re (F k N k ) > Id . 

4. Choice of F^. This is the new part detailed in the Appendix A of 
|MZj . Denote by £ k the block diagonal matrix T)i&g(E k ). A vector w E C^ fe , 
N k = v k a kl is broken into a k blocks w p £ C^ fc , with components denoted by 
w P) a- We denote by the v k x v k blocks of lZ k and by R Pt a,q,b their entries. 
The entries of E k are denoted by E a ^. Since R pA ,q,b = when b > 1, the 
special form of E k implies that 

Re(£ k 1l k W,w) = Re 'Y^ E a,cRp,a,q,lWq,XWp,c 

= Re ^ ek,ir P ,qWq,iw p ,i + 0(\w^i\ |u>'J) 



19 



where G C ak is the collection of the first components w Pt i, w'^ the 

remainder components and r P)IJ = R PtUk ,q,i the lower left hand corner entry 
of R p ,q- The matrix Re(R k ) is definite, positive or negative according to 
the sign of q k , which is the sign of e k ,i by (|5.'29|) . Thus, multiplying E k by 
some positive constant, we can achieve that in addition to (|5.30|) (|5.32|) . the 
following inequality holds: 

(5.33) Re(S k Kk(pX,0)w,w) > 2\w^\ 2 - C'\w'J 2 . 

for some constant C > 0. 

As in 3, there is Fl real and skew symmetric such that for all w € C fe : 

Re(F^N k w,w) > -\w x \ 2 + {C + l)\w'\ 2 . 

Thus, with T' k = Diag(i^), M k = Diag(Af fc ) and weC N <°: 

Re{F' k N k w,w) > -\w* A \ 2 + (C + l)\w'J 2 . 

Therefore, with (|5,33|) . we have 

(5.34) Re (£ k K k (p,(,0) - iT k Q k (jp,Q) > Id. 

5. Summing up, we see that the matrix defined in (|5.27f) is self adjoint. 
From (|5.30j) . the condition (|5.14j) is satisfied provided that 5c = 1. Moreover, 

Re {S k {Q k + P K k ))= ^8V k (p, C) + p5V' k ( P , (, p), 

with T> k = Diag(Dfc) and, at the base point, 

V' k (p,(,0) = Re (£ k n k (p,C,0) - iT' k Q k (p,()) . 

By 1)5.32(1 and 1)5.34(1 . the matrices V k and V k are positive definite on a 
neighborhood of the base point (p, (, 0). This implies that S k and 7i\ satisfy 
(t5T5|) . 

This finishes the construction of symmetrizers for TL\ and thus the proof 
of Theorem 14.11 
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